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1.  INTRODUCTION 


* 


The  introduction  of  the  Fluidyne  “no- moving- parts”  pump  in  197 1 12  stimulated  some  new 
interest  in  the  conversion  of  low  grade,  low  temperature  heat  energy  into  useful  work.  The  basic 
principle  was  reported  to  be  a closed  Stirling  cycle,  wherein  two  adjacent  U-tubes  of  liquid  were 
stimulated  into  oscillation  by  the  application  of  heat  to  a connecting  U-tube  of  gas.  The  open 
U-tube  of  oscillating  liquid  thus  was  the  source  of  power  to  do  work.  A linearized  lossless  model 
was  proposed  by  Elrod’  that  predicted  the  onset  of  oscillation  and  the  frequency  thereof.  No 
satisfactory  physical  explanation  was  offered  for  the  phenomenon,  however. 

With  the  coming  of  the  energy  crisis,  initial  studies  on  reclaiming  waste  heat  were  proposed 
by  the  Harry  Diamond  Laboratories  using  the  Fluidyne.  With  the  proposed  increased  dependence 
on  nuclear  power,  additional  interest  in  the  Fluidyne  became  apparent  when  it  was  conceptually 
demonstrated  that  a Fluidyne  could  serve  as  an  independent  power  source  for  fluidic  monitoring 
and  control  schemes  that  are  completely  decoupled  from  standard  electromechanical  systems  in 
use.  Such  incidents  as  the  Brown’s  Ferry  fire  and  the  loss  of  secondary  controls  could  be  avoided  by 
using  the  residual  nuclear  heat  to  operate  a self-contained  Fluidyne  pump  and  fluidic  sensors  and 
controls  operating  within  the  pressure  vessel  itself.  A system  having  no  moving  parts  except  the 
working  fluid  would  be  extremely  reliable. 

The  poor  efficiency  of  currently  available’  pumps  allows  profitable  operation  only  where  the 
heat  source  costs  are  negligible  or,  as  in  a nuclear  plant,  not  of  concern.  Improvement  of  the  effi- 
ciency would  allow  greater  utilization  and  would  decrease  the  unit  size  per  unit  output  power. 
This  decrease  would  allow  the  use  of  smaller  devices  in  such  applications  as  nuclear  power  plant 
auxiliary  control  systems. 

This  report  presents  a comprehensive  basis  for  the  understanding  of  the  Fluidyne  pump’s 
operation,  analyzes  the  concept  of  the  efficiency  of  the  thermodynamic  cycle,  and  provides  impor- 
tant guidelines  for  design. 


2.  PHYSICAL  DESCRIPTION  OF  FLUIDYNE  STARTUP 

In  the  schematic  diagram  (fig.  1),  a reservoir  of  liquid  is  shown  with  a U-tube  of  gas  con- 
necting two  liquid  columns.  A separate  U-tube  filled  with  the  liquid  is  situated  with  one  end 
directly  beneath  one  end  of  the  gas-filled  tube.  The  other  end,  for  this  analysis,  is  left  open  to  the 
atmosphere.  The  free  surface  motion  indicates  the  output  mechanical  work  available. 

The  gas  volume  may  be  thought  to  be  composed  of  two  imaginary  halves.  The  right  half, 
labeled  “2,”  which  connects  to  the  output  tube,  is  the  hot  side;  the  left  half,  labeled  “1,”  is  the  cold 
side.  Heat  is  continually  applied  to  the  right  half,  and  the  left  side  is  maintained  at  a lower  or 
ambient  temperature.  In  the  experiments  conducted,  the  cold  side  was  cooled  by  natural  convec- 
tion to  ambient. 

■C.  West,  The  Fluidyne  Heat  Engine,  United  Kingdom  Atomic  Energy  Authority.  Harwell,  UK,  AERE-R-6775  (1971). 

lC.  West,  And  Yet  It  Moves,  New  Scientist,  29  (August  1974). 

‘H.  Elrod.  The  Fluidyne  Heat  Engine-How  to  Build  One,  How  it  Works,  National  Technical  Information  Service. 
Springfield,  VA,  AD  A006367  (1974). 


5 


mm 


MOTION  OF  UOUIO 


LIQUID  RESERVOIR 


Figure  1.  Schematic  of  Fluidyne  pump. 
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Figure  2.  Schematic  representation  of 
Fluidyne  startup  process. 


a.  In  equilibrium  conditions,  no  heat  is  applied.  The  liquid  levels  h,,  h,,  and  h»  are 
prescribed  by  the  static  pressure  in  the  gas  volume  of  columns  (1)  and  (2),  and  h,  = h,.  (If  the  gas 
pressure  p,,  is  ambient,  p,,,,,  then  all  three  levels  are  equal.) 


! 
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b.  Heat  is  applied  to  the  right  half  of  the  gas  volume  of  column  (2).  The  gas  expands,  and 
the  gas  volume,  V reases  while  the  gas  pressure  p,2  also  increases.  As  a natural  consequence, 
liquid  levels  h,  and  ecrease,  and  liquid  level  h5  increases  as  the  d in  the  displacer  is 
accelerated.  The  gas  temperature  rises  in  column  (2). 

c.  The  inertia  of  the  liquid  in  the  output  tube  keeps  the  liquid  moving;  hence,  the  liquid  in 
column  (2)  tends  to  keep  moving  downward,  also.  As  more  hot  side  is  exposed,  the  liquid  is 
accelerated  further  until  equilibrium  in  the  gas  is  reached. 

d.  Gravity  stops  the  upward  motion  of  liquid  level  hi,  and  levels  h2  and  hs  tend  to  seek 
equality.  Hence  h2  increases,  and  h5  starts  to  decrease.  At  the  same  time,  there  is  a reduction  of  hot 
volume.  This  reduction  starts  an  overall  gas  volume  decrease  due  to  a temperature  reduction  since 
less  gas  is  exposed  to  the  heat  source.  The  pressure  p,2  decreases.  All  this  action  is  augmented  by  the 
descent  of  the  output  level  due  to  gravity,  which  produces  a dynamic  head  at  the  base  of  column 
(2),  further  increasing  h2.  An  overshoot  of  the  two  liquid  levels  occurs,  h2  > ht,  h2  > h9.  This  over- 
shoot gives  rise  again  to  a gravity  potential  between  liquid  levels. 

e.  The  liquid  levels  again  seek  equilibrium  with  level  h2  decreasing.  This  decrease  exposes 
more  of  the  gas  to  the  heat  source  and  causes  the  gas  to  expand;  the  cycle  begins  all  over  again, 
with  the  output  liquid  again  moving  upward. 

This  is  the  physical  description  as  observed  in  a glass  or  plastic  demonstration  model  built  and 
tested  during  this  program. 


r 


t 
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3.  MATHEMATICAL  MODEL 

Geisow4  has  proposed  a set  of  governing  equations  for  a lossless  Fluidyne  system  in  an  attempt 
to  compute  the  conditions  required  for  the  onset  of  instability.  To  this  end,  he  used  the  Euler  equa- 
tions of  motion  in  one  dimension,  linearized  them,  and  then  performed  a first-order  perturbation 
stability  analysis.  The  results  of  this  analysis  led  to  an  expression  for  the  critical  temperature  dif- 
ferential, AT,  as  a function  of  the  pump  geometry  and  for  the  oscillating  frequency,  oj.  This  is 
given  as  equations  (1)  and  (la): 


(1) 


• (la) 


'A.  D.  Geisow.  The  Onset  o)  Oscillations  in  a Lossless  Fluidyne,  United  Kingdom  Atomic  Energy  Authority,  Harwell, 
UK.  AERE-M-2840  (October  1976). 
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where 


Tj  = cold  temperature  (K) , 

T2  = hot  temperature  (K) , 

p = liquid  density  (kg/m3), 

g = acceleration  of  gravity  (m/s2) , 

= equilibrium  volume  of  gas  (m3) , 

= cross-sectional  area  of  displacer  tube  (m2) , 

= atmospheric  pressure  (Pa) , 

£ = equilibrium  length  of  cold  leg  to  output  (m) , 

c 

= equilibrium  length  of  hot  leg  to  output  (m) , 

= length  of  displacer  (m), 

Ac  - cross-sectional  area  of  cold  leg. 


This  analysis  was  performed  for  the  geometry  shown  in  figure  3 and  indicates  a strong 
dependency  of  system  response  on  the  displacer  tube. 


Figure  3.  Schematic  of  Fluidyne  II  used  by  Geisow  (AERE-M-2840, 
October  1976). 


The  original  Fluidyne  proposed  by  West'  and  described  by  Elrod3  did  indeed  have  such  a con- 
figuration; however,  currently  commercially  available  Fluidynes5  replace  the  displacer  tube  with 
a fluid  reservoir  as  shown  in  figures  1 and  4. 

'C.  West,  The  Fluidyne  Heat  Engine,  United  Kingdom  Atomic  Energy  Authority.  Harwell,  UK,  AERE-R-6775  (1971). 

'H.  Elrod,  The  Fluidyne  Heat  Engine-How  to  Build  One,  How  it  Works,  National  Technical  Information  Service, 
Springfield,  VA,  AD  A006367  (1974). 

'Stirling-Cycle  Engines  Developed  at  Harwell-The  Fluidyne  3 Pump.  United  Kingdom  'Atomic  Energy  Authority. 
Harwell.  UK.  Harwell  Design  Studios  56/1975  (1975/ 
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Figure  4.  Schematic  of  actual  Fluidyne  pump  (diode  elements— ball 
valves — permit  unidirectional  flow). 


Interestingly  enough,  the  models  shown  by  Elrod3  always  have  the  output  tube  axially  in  line 
with  the  hot  leg  as  opposed  to  being  at  right  angles  as  shown  in  figure  3,  even  when  the  displacer  is 
present.  This  leads  to  the  observation  of  section  2 that  dynamic  head  feedback  from  the  output  to 
the  hot  leg  is  important  in  the  pump  operation  in  that  the  inertia  of  the  moving  fluid  in  the  output 
leg  preferentially  acts  to  increase  the  head  in  the  hot  leg. 

The  following  set  of  governing  equations  may  be  generated,  assuming  slug  flow  in  all  tubes 
for  the  geometry  of  figure  1 or  4. 

The  pressure  drop,  APjj  =Pj  - Pj,  for  any  tube  is  given  by  the  electrical  analog  equation  for 
time  varying  slug  flow*  (derived  from  the  Euler  equation): 


AP . . 

ID 


dQ.  . 

L - 
ij  dt 


R. .Q.  . + h.p  g 
ID  ID  i L 


(2) 


where 


L.  . 

ID 


ID 


ID 


liquid  inductance  (kg/m4)  between  points  i and  j,  = j 

liquid  density  (kg/m3) , 

tube  length  between  i and  j (m) , 

tube  cross-sectional  area  (m'), 


>H.  Elrod.  The  Fluidyne  Heat  Engine-How  io  Build  One,  How  it  Works,  National  Technical  Information  Service, 
Springfield.  VA.  AD  A006367  (1974). 

*J.  M.  Kirshner , Fluid  Amplifiers,  McGraw-Hill  Book  Co,  New  York  (196 7). 
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®i j " volumetric  flow  rate  (m3/s)  in  direction  i to  j, 
t = time  (s) , 

Rij  = li(luid  resistance  (for  laminar  flow  R = 8Ttuf/A2)  [kg/tm4  • s)]  , 
M = liquid  viscosity  [kg/(m-s)] 
hi;j  = liquid  height  (m)  , 
g = acceleration  of  gravity  (m/s2) . 


Therefore,  for  the  cold  leg,  between  steps  (a)  and  (c)  the  pressure  drop  is 


p3  ~ P]  l13  dt  + R31Q31  + hlPig  . (3) 

Similarly,  the  hot  leg  pressure  drop  can  be  written  when  the  additional  total  pressure  feedback 
term  is  included  due  to  the  dynamic  head  effect  of  the  output  leg: 


P4  ' P2  - L4  2 + R4  2®42  + h2PLg 


■MS'  ■ 


where 


K1  = °'  ®4  5 > °»  Kl  = 1.  £>45  < 0. 


Finally,  the  output  leg  presure  drop  is,  with  hot  leg  feedback. 


P4  - P5  = 


dg4  5 PL  /Q42\2 

dt  + R45  245  + h5PLg  - — K2^— J , 


where 


K2  = 0,  Q42  > 0;  K2  = 1,  g42  < 0. 


. . . :i 


The  equation  relating  the  liquid  height  and  the  volume  flow  rate  is  the 


constitutive  equation 


a dh  ^ 

A dF  ’ C • 


(6) 


Hence,  the  individual  time  rates  of  height  change  are 


dt 

dh2 

dt 

dh 


‘31 
^3  1 


*4  2 


42 


5 _ 2i*5 

^ " \7 


(7) 


(8) 


(9) 


It  follows,  from  conservation  of  mass  of  an  incompressible  liquid  into  the  reservoir,  that  the  net 
tlow  must  be  zero,  or 


- Q - Q -O  =0 
31  *42  *45 


(10) 


1 


Note  the  minus  signs  since  each  flow  is  assumed  to  be  positive  out  of  the  reservoir.  The  pressure  at 
the  base  of  the  cold  leg  is  the  reservoir  pressure.  If  the  output  leg  does  not  impede  flow  from  the 
reservoir  into  it  or  the  hot  leg,  then  pj  = p4.  If  there  is  an  orifice  resistance,  then 


PL«31 


3 " 2(A42  - A45)2  V ' 

where  cj  is  the  discharge  coefficient. 

For  convenience  of  this  initial  analysis,  consider  p,  to  be  ambient: 

P5  = 101  kPa  (absolute) 


(ID 


(12) 


We  have  up  to  this  point  9 equations  and  11  unknowns— namely,  5 pressures,  3 flows,  and  3 
heights;  consequently,  more  equations  are  required. 

If  the  assumption  is  made  that  the  pressure  is  constant  at  any  instant  of  time  in  the  gas 
volume,  which  is  a reasonable  assumption  considering  that  the  process  is  relatively  slow,  the 
gaseous  resistance  of  the  tubing  is  negligible,  and  the  head  potential  is  likewise  negligible,  then 
another  equation  is  available: 


P 


l 


(13) 


The  gas  mass  must  be  conserved  so  that 


nij  + m2  = constant  . (14) 

This  equation  has  added  two  unknowns,  mi  and  m2,  the  masses  of  the  cold  and  hot  sides  of  the 
gas  volumes  (kg) , respectively;  hence,  there  are  now  1 1 equations  and  13  unknowns.  The  equation 
of  state  may  be  written  for  each  side  of  the  gas  volume.  If  it  is  assumed  that  the  gas  is  ideal,  then 


PlV! 

= mjiTTj  , 

(15) 

PlV2 

= m2i?T2  , 

(16) 

where 

Vj  = volume  of  cold  gas  (m3) , 

V2  = volume  of  hot  gas  (m3) , 

Tj  = temperature  of  cold  gas  (K) , 

= temperature  of  hot  gas  (K) , 

/?  = gas  constant  [m3/(s2K)J 

There  are  now  13  equations  and  17  unknowns,  the  new  unknowns  being  volumes  V,  and  V2 
and  temperatures  Tt  and  T2. 

The  gas  volumes  are  directly  related  geometrically  to  the  position  of  the  gas-liquid  interface  so 

that 


'l  =K  ' (h!  - he)A13 


(17) 


12 


and 


'2  -K  - (*2  - »,)\ 


where 


hg  = equilibrium  height  of  liquid  (m) 


This  now  brings  us  to  15  equations  and  17  unknowns. 

At  this  point,  a statement  can  be  made  about  temperatures  T,  and  T2.  Elrod3  and  GeisoW 
chose  to  assume  that  the  temperatures  are  known  constants.  This  choice  is  tenuous  since  it  requires 
that  a step  distribution  in  temperature  exists  in  the  gas  tube.  For  the  sake  of  simplicity,  assume  that 


Tj  = constant  , 
T2  = constant  , 


where  T,  and  T2  are  average  temperatures  (K)  in  each  volume  of  gas.  These  then  bring  the  count  of 
unknowns  and  equations  even  at  17  apiece.  The  evaluation  of  T,  and  T2  is  given  in  appendix  A. 
Upon  examination  of  the  governing  equations,  one  notes  that  the  pressure  drop  equations  are 
nonlinear  since  the  inductance  and  the  resistance  terms  are  functions  of  the  length  of  the  liquid, 
which  is  variable.  Linearization  can  be  achieved  by  stipulating  that  the  total  excursions  of  the 
liquid  interfaces  be  considerably  smaller  than  the  equilibrium  lengths.  For  computer  solution, 
however,  this  stipulation  is  not  necessary.  There  are  five  independent  first-order  ordinary  differen- 
tial equations,  since  one  of  the  flow  equations  is  used  to  find  a pressure.  The  solution  of  the  entire 
set  of  equations  gives  the  displacement  of  the  output  column  of  liquid.  Since  the  object  of  this  study 
is  to  establish  the  parameters  relating  to  the  efficiency  of  the  pump,  a quantity  of  interest  is  the 
output  power,  W.  This  is  the  time  rate  of  increase  of  potential  energy  of  the  output  or  the  root 
mean  square  (rms)  value  of  the  product  of  the  flow  and  the  dynamic  head.  The  instantaneous 
power  out  is 


SH.  Elrod,  The  Fluidyne  Heat  Engine-How  to  Build  One,  How  it  Works.  National  Technical  Information  Service 
Springfield,  VA,  AD  A006367  (1974). 

'A.  D.  Ceisow,  The  Onset  of  Oscillations  in  a Lossless  Fluidyne,  United  Kingdom  Atomic  Energy  Authority.  Harwell 
UK,  AERE-M-2840  (October  1976). 


The  useful  power  is  the  rms  value  or 


W = 0.707W  . (22) 

rms  max 

The  input  energy  is  the  heat  input  from  which  the  temperature  distribution  was  obtained 
(app  A). 

Physical  observations  indicate  that  most  of  the  losses  occur  in  the  input  of  heat  to  the  system. 

Equations  (3)  to  (20)  may  now  be  cast  into  a computer  program  and  solved  simultaneously 
without  resorting  to  perturbation  or  threshold  analyses.  This  simultaneity  is  important  because  the 
actual  values  of  displacements  are  needed  to  compute  the  output  power. 

4.  NUMERICAL  RESULTS 

The  computer  program  is  given  in  appendix  B with  a sample  printout.  Figures  5 to  8 plot  the 
various  parameters  of  interest  for  typical  geometries  in  which  oscillations  are  self-sustained. 

In  figure  5,  the  three  liquid  displacements,  the  gas  pressure,  and  the  change  of  gas  volume  are 
shown  as  a function  of  time.  They  are  in  substantial  agreement  with  the  descriptive  model 


TIME  III 

Figure  5.  Time  history  of  Fluidyne  response:  L„  = liquid  inductance  output,  Ve  = volume 
equilibrium,  A„  = output  area,  Tj  = temperature  of  hot  side  of  gas-liquid  interface,  and  h = 
height  of  liquid. 


proposed  as  a result  of  physical  observations.  Constant  temperature  apparently  causes  the  oscilla- 
tions to  keep  growing.  For  limit  cycling,  a more  comprehensive,  flow  dependent  temperature 
distribution  has  to  be  developed  that  in  effect  provides  a time  delay  in  the  temperature  change. 
The  heat  conduction  time  constant  for  air  is  t = £>CprVk  = 7 s,  where  Cp  = specific  heat  at  constant 
pressure  (W  • s/kgK) , r = radius  of  the  pipe  (m) , and  k = thermal  conductivity  (W/mK) . Hence,  it  is 
a poor  assumption  that  the  temperature  reaches  the  hot  or  cold  value  instantaneously,  especially 
since  the  period  of  oscillations  is  an  order  of  magnitude  less  than  7.0  s.  Similarly,  the  convective 
time  constant  due  to  the  movement  of  the  gas-liquid  interface  must  be  of  the  same  order,  that  is, 
1 s.  Due  to  the  complexity  of  the  fluid  mechanic-heat  transfer  governing  equations  for  such  a 
phenomenon,  this  development  must  be  left  to  a more  comprehensive  and  detailed  analysis  at  a 
later  time. 


In  figure  6,  one  sees  the  variation  of  amplitude  and  frequency  of  oscillation  for  different 
parametric  changes.  If  one  observes  the  lower  three  sets  of  curves,  the  effect  of  increasing  the  out- 
put length  can  be  seen.  As  one  would  expect,  the  period  increases  with  output  length  due  to  increased  inductance. 


Figure  6.  Output  displacements  for  various  parameters:  LQ  = liquid  inert- 
ance  output,  T2  - temperature  of  hot  side  of  gas-liquid  interface,  Ve  = 
volume  equilibrium,  and  AQ  = output  area. 
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The  middle  set  of  curves  demonstrates  the  effect  of  increasing  the  hot  side  temperature.  The 
frequency  decreases  slightly  with  increasing  temperature.  This  decrease  is  a result  of  the  higher 
mean  pressure  in  the  gas  volume  that  results  in  a larger  volume  and,  hence,  a slightly  increased  gas 
capacitance.  Interestingly,  the  amplitude  appears  to  have  a maximum  at  a midrange  temperature. 

The  curve  marked  “305”  shows  the  effect  of  increasing  the  output  area.  This  increase  has  the 
effect  of  decreasing  the  inductance,  again  increasing  the  frequency. 

The  top  curve  of  figure  6 shows  the  effect  of  increasing  the  inital  gas  volume  and  the  mass. 
The  gas  volume,  as  already  suspected,  acts  as  a capacitance;  hence,  increasing  the  capacitance 
decreases  the  frequency,  as  demonstrated  in  the  figure. 

Intuitively,  one  would  expect  that  the  predominant  frequency  determining  terms  are  the  out- 
put inductance  and  the  gas  capacitance.  This  would  lead  to  the  assumption  that  the  system  fre- 
quency, f,  should  behave  approximately  as  an  RLC  circuit,  so  that 


f 


2,(Vg)^ 


(23) 


where 


CG  = gas  capacitance  = ve^P„  (m4s2/kg)  , 

p^  = equilibrium  absolute  gas  pressure  (Pa)  . 

Table  I shows  the  comparison  of  the  frequencies  from  equation  (23)  with  the  actual  numerical 
results  obtained  from  the  curves  in  figure  6. 


TABLE  I.  COMPARISON  OF  NUMERICALLY  COMPUTED  FLUIDYNE 
FREQUENCY  WITH  EQUATION  (23) 


Parameter 

Frequency 

(Hz) 

Output 

Gas 

Output 

Case 

length 

volume 

area 

Numerical 

Eq  (23) 

No. 

(m) 

(xlO-5  m3) 

(xlO'4  m3) 

2.0 

2.0 

1.0 

2.50 

2.5299 

1 

7.0 

2.0 

1.0 

1.54 

1.3520 

2 

20.0 

2.0 

1.0 

0.77 

0.7997 

3 

2.0 

4.0 

1.0 

1.74 

1.7883 

4 

2.0 

2.0 

2.0 

3.33 

3.5766 

5 

Note:  Liquid  density  = 1000  kg/m1;  ambient  absolute  pressure  p^  - 101  kPa. 


r * 


■ » 


The  agreement  is  good,  so  equation  (23)  may  be  used  as  a design  guide  for  determining  the 
operating  frequency  of  a Fluidyne. 

The  question  of  cycle  efficiency  can  be  answered  only  by  determining  the  nature  of  the  ther- 
modynamic cycle,  the  energy  input,  and  the  energy  output.  In  figure  7,  the  relation  between  gas 
pressure  and  volume  is  shown  for  startup  and  the  first  cycle  for  case  1 of  table  I.  It  is  immediately 
apparent  that  this  is  not  a Stirling  cycle  as  originally  thought  by  the  inventor.'  Instead,  it 


2 kPa  x 0.5x10~*m3  . kPa_m3  „ jqj  P« 
m2  m2 


« 10  3 


Pam3 


pV  * 0 525  x 10  3P*  m3 
X-  1/1 -0.4* 

’3,“«’03w 

RMS  OUTPUT  POWER  FIRST  CYCLE  - (°  ^ * ° U2  ) * 1fr3W  * 0 707 


12.36% 


- 0 16225  * 10  3W 


resembles  more  a Carnot  cycle,  albeit  it  is  not  that,  either.  The  observed  efficiency  of  the 
Fluidyne’s  first  cycle  is  12.36  percent  when  defined  as  the  output  power  divided  by  the  power  of 
one  cycle  (/ p dV) . The  heat  input  to  the  gas  gives  the  overall  efficiency.  The  heat  input  may  be  ap- 
proximated by  the  convective  heat  transfer  equation  and  the  mean  gas  temperature  Tj,  as  com- 
puted in  appendix  A: 


'C.  West,  The  Fluidyne  Heat  Engine,  United  Kingdom  Atomic  Energy  Authority,  Harwell,  UK,  AERE-R-6775  (1971). 
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(24) 


where 


H = convective  heat  transfer  coefficient  (W/K/m2) , 

A = outside  area  of  hot  volume  = 2irri  (m2)  , 

G 

T = applied  hot  temperature  (K) 
h 

For  the  case  in  question,  the  average  gas  temperature  Tj,  for  an  applied  temperature  T},  = 
315  K is  313.454  K.  An  average  value  of  H = 17  W/K/m2,  and  the  gas  hot  side  area  is  3.77  * 10~3  m2. 
These  values  result  in  a heat  input  of  0.21  W and,  consequently,  a poor  efficiency  overall.  When 
the  pump  is  allowed  to  operate  for  three  cycles,  whereupon  it  is  expected  that  almost  steady-state 
conditions  would  normally  have  been  reached,  the  overall  efficiency  is  up  to  1.2  percent.  The 
third-cycle  efficiency  is  shown  as  a function  of  the  applied  differential  temperature  in  figure  8.  For 
L„  = 2.0  m,  there  is  a maximum  at  a differential  temperature  of  about  24  K.  This  maximum  could 
have  been  predicted  by  the  observations  from  figure  6,  where  it  was  shown  that  the  amplitude  of 
oscillation  was  lower  at  elevated  temperatures. 


■77  a RMS  OUTPUT  POWER  = 0.707W nu» 
HEAT  INPUT  hA^T 


Figure  8.  Fluidyne  efficiency  as  function  of  input 
temperature  and  output  length:  h = height  of  liquid 
column  and  L0  = liquid  inertance  output. 
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Also,  it  can  be  observed  in  figure  6 that  the  amplitude  of  oscillation  was  inversely  propor- 
tional to  the  output  length;  therefore,  it  seems  intuitive  that  the  efficiency  increases  if  the  output 
length  decreases.  When  the  gas  volume  increases,  the  amplitude  also  increases,  but  the  frequency 
decreases.  The  efficiency  may  increase;  however,  since  the  frequency  is  lower,  the  output  flow  is 
lower,  so  the  net  effect  is  that  the  efficiency  is  not  materially  affected.  Computed  results  (not 
shown  here)  bear  out  this  thesis.  Figure  8 shows,  however,  that  efficiency  is  markedly  improved  by 
decreasing  the  output  length.  Also  shown  is  that  the  optimum  differential  temperature  is  lowered 
by  a factor  of  two  as  the  maximum  efficiency  is  more  than  doubled. 


5.  DISCUSSION  AND  CONCLUSIONS 

The  disturbing  point  in  this  analysis  is  the  exponential  growth  of  oscillations.  Appendix  A 
discusses  this  growth  and  the  temperature  distribution.  More  fundamentally,  however,  the  reason 
for  this  growth  is  probably  the  assumption  that  the  gas  changes  temperature  instantaneously  from 
hot  to  cold  or  vice  versa.  It  takes  an  infinite  amount  of  power  or  energy  to  raise  a finite  amount  of 
mass  a finite  temperature  in  zero  time.  In  a numerical  model,  this  infinity  is  tempered,  but  still  re- 
mains large.  The  conclusion  to  be  drawn,  therefore,  is  that  the  gas  must  be  much  more  carefully 
modeled. 

The  values  of  efficiency  presented  here  cannot  be  used  in  anything  more  than  a comparative 
sense  since  they  are  based  on  output  power  values  that  are  steadily  increasing.  Observations, 
however,  have  indicated  that  the  actual  process  reaches  a steady  state  after  about  three  cycles. 
Since  the  frequencies  computed  are  close  to  those  observed  and  the  amplitudes  are  reasonable,  the 
three-cycle  numerical  results  have  been  used. 

Two  efficiencies  are  involved,  thermodynamic  and  mechanical.  Thermodynamic  efficiency  is 
relatively  low,  and  mechanical  efficiency,  although  an  order  of  magnitude  greater  than  thermo- 
dynamic, is  still  only  about  10  percent.  Thus,  the  low  efficiencies  reported  are  substantiated  by 
this  numerical  analysis.  There  is  hope,  however,  of  improving  the  overall  efficiency  by  a factor  of 
two  or  even  three  over  the  extant  devices  by  a thorough  parametric  study.  This  study  has  examined 
more  than  2500  permutations  of  geometry,  and  they  cannot  be  presented  here.  The  important 
observations,  however,  have  been  presented,  and  hope  for  better  efficiency  is  offered. 

Several  conclusions  may  be  drawn  from  this  study.  The  Fluidyne  pump  operates  not,  as  first 
thought,  on  a Stirling  cycle,  but  on  a cycle  more  closely  resembling  a Carnot  cycle.  It  is  essentially 
a low-efficiency  device,  although  overall  thermal  efficiencies  of  10  percent  may  conceivably  be 
reached.  Poor  thermal  coupling  with  the  gas  accounts  for  a large  portion  of  this  inefficiency.  In  the 
real  device,  it  would  be  advantageous  to  have  a nonconducting  liquid  and  a highly  conducting  gas 
(a  contradiction  of  thermodynamic  terms) . This  advantage  could  possibly  be  realized  by  providing 
an  insulating  interface  substance  such  as  foam.  The  gas  volume  should  be  insulated  so  that  all  the 
heat  input  stays  in  the  gas. 

Relating  to  the  analysis  itself,  the  important  points  and  parameters  have  been  identified.  The 
displacer  of  the  original  Fluidyne  has  been  shown  to  be  superfluous  to  the  device.  Although 
amplitude  oscillations  were  never  constant,  limit  cycling  might  occur  when  the  temperature 
distribution  is  corrected  or  modeled  more  appropriately.  A simple  model  for  the  frequency  of  the 
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Fluidyne  has  been  proposed  that  relates  the  gas  volume  to  the  geometry  of  the  output  tube  in  terms 
of  an  electrical  RLC  second-order  system  analogy.  Good  agreement  with  numerical  results  has 
been  demonstrated. 

The  overall  objectives  of  this  study  were  to  determine  if  the  efficiency  of  existing  Fluidynes 
could  be  increased  and  to  comprehensively  describe  the  physics  of  the  Fluidyne  operation.  These 
objectives  have  been  met.  The  efficiency  may  be  increased  by  appropriate  geometric  choices,  and 
for  the  first  time  a clear  description  has  been  proposed  for  the  operation. 

Research  and  development  of  the  Fluidyne  heat  engine  should  continue.  The  Fluidyne  can 
provide  fluid  power  with  tiny  temperature  differentials,  with  no  moving  parts.  For  power  produc- 
tion in  inaccessible  locations,  this  device  is  ideal.  One  of  the  most  inaccessible  places  is  the  interior 
core  of  the  pressure  vessel  of  a nuclear  reactor.  The  nuclear  community  requires  redundant  and 
completely  independent  safety,  monitoring,  and  control  systems.  Presently,  all  backup  controls  de- 
pend on  a single  type  of  energy  source — electricity.  Should  this  ever  fail,  then  no  systems  other 
than  the  purely  mechanical  ones  could  operate.  The  mechanical  systems  normally  can  provide  on- 
ly shutdown  capability.  A fluidic  system,  self-contained  within  the  reactor  core,  with  no  pressure 
vessel  penetrations,  powered  by  the  rectified  output  of  a Fluidyne  can  provide  all  the  active 
monitoring  and  control  of  the  standard  systems  and  also  cause  shutdown  in  the  event  of  a scram. 
This  approach  can  avoid  the  costly  and  inconvenient  automatic  shutdowns  presently  envisioned.  A 
reactor  may  continue  to  operate  with  the  fluidic  system  installed. 

Recent  panel  discussions  sponsored  by  the  American  Society  of  Mechanical  Engineers 
( ASME)  * on  fluidics  and  nuclear  engineering  specifically  mentioned  fluidic  shutdown  schemes  be- 
ing developed  in  the  United  States.  At  the  same  session,  fluidic  two-diode  pumping  systems  were 
described  that  have  been  developed  for  the  United  Kingdom  Atomic  Energy  Authority.  These 
fluidic  diode  systems  have  successfully  operated  with  Fluidyne  pumps  at  Harwell.  C.  F.  King 
of  the  University  College,  Cardiff,  UK,  reported  that  newly  developed  vortex  diodes  that  have 
turndown  ratios  in  excess  of  30:1  operate  exceptionally  well  when  attached  to  the  Fluidynes.  For  a 
session  on  energetics,  fluidic  temperature  sensors  that  can  withstand  extremely  hot,  harsh  en- 
vironments were  reported  by  Drzewiecki  et  al.7  Such  sensors  are  ideal  candidates  for  use  with  a 
Fluidyne  system  for  the  monitoring  of  nuclear,  fossil  fuel,  solar,  and  geothermal  power  generation 
plants.  Of  particular  interest  are  military  nuclear  power  facilities.  These  are  normally  mobile  and 
do  not  have  a ready  source  of  auxiliary  power  to  operate  secondary  systems.  Even  when  they  do, 
the  proximity  of  such  auxiliary  systems  makes  them  vulnerable  to  failure  due  to  either  the  military 
environment  or  the  radiation  possibility  from  the  primary  plant  itself.  The  Fluidyne  offers  a 
unique  solution  to  this  problem.  It  merely  becomes  an  integral  part  of  the  whole. 


'T  M.  Drzewiecki,  Fluidic  Applications  in  Nuclear  Engineering,  ASME  Winter  Annual  Meeting,  Atlanta,  GA 
(November  1977). 

’T.  M Drzewiecki,  R.  M.  Phillippi.  and  C.  Paras,  Fluidics-A  New  Potential  for  Energy  Conservation  by  Continuous 
High  Temperature  Monitoring  and  Control,  Heat  Transfer  in  Energy  Conservation,  R.  J.  Goldstein  et  al.  ed.,  ASME  Book 
H00106.  New  York  (November  1977),  127-134. 
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Other  military  applications  for  remotely  powered  fluidic  systems  are  battlefield  security  and 
sensor  systems  (REMBASS).  Fluidic  sensors  have  been  demonstrated  to  be  uniquely  suited  to  the 
battlefield.  The  source  of  power  has  been  of  concern,  however.  A Fluidyne  that  can  operate  on  a 
3-  or  4-deg  temperature  difference  could  operate  during  the  day  on  light  energy  and  during  the 
night  on  the  stored  heat  in  the  ground  with  the  colder  night  air  being  the  thermal  sink.  This  discus- 
sion is  limited  only  by  imagination.  It  is  the  author’s  opinion  that  many  military  systems  problems 
can  be  effectively  and  efficiently  solved  by  fluidic,  flueric,  or  fluid  systems  powered  by  Fluidyne 
pumps. 
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NOMENCLATURE 


A area  (m2) 

c discharge  coefficient  (-) 
d 

c specific  heat  at  constant  pressure  (W-s/kgK) 

P 4 - 

C capacitance  (m  sVkg) 
f frequency  (Hz) 
g acceleration  of  gravity  (m/s2) 
h height  of  liquid  column  (m) 

H convective  heat  transfer  coefficient  (W/K/m2) 
k thermal  conductivity  (W/mK) 

K constant,  =1  when  Q < 0;  =0  when  Q > 0 
l column  length  (m) 

L liquid  inertance  (inductance)  (kg/m1*) 
m mass  (kg) 

p absolute  pressure  (Pa) 

P pressure  (Pa) 

Q volumetric  flow  rate  (m3/s) 
r radius  of  pipe  (m) 

R liquid  resistance  [kg/(m4-s)] 

R gas  constant  [m3/(s2K)j 
t time  (s) 

T absolute  temperature  (K) 

T mean  gas  temperature 
u liquid  velocity  (m/s) 

V volume  (nr ) 

W power  (W) 

A change 

p fluid  absolute  viscosity  [kg/(m-s)] 
n efficiency 
p liquid  density  (kg/m3) 

T time  constant  (s) 

u radian  oscillating  frequency  (rad/s) 


APPENDIX  A 


The  gas  volume  that  provides  the  means  for  the  motive  power  in  the  Fluidyne  is  essentailly 
cylindrical.  It  is  bounded  on  either  end  by  liquid  interfaces  that  move  so  that  the  volume  of  gas  ex- 
pands and  contracts.  Heat  is  applied  to  one  half  and  rejected  from  the  other  half  by  maintaining 
essentially  constant  hot  and  cold  wall  temperatures,  Th  and  Tc,  respectively.  Consider,  therefore, 
two  touching  cylinders  with  constant  wall  and  end  temperatures.  The  end  temperatures  can  be 
justified  to  be  Th  or  Tc  if  the  conductivity  of  the  liquid  is  much  greater  than  that  of  the  gas.  Figure 
A-l  shows  the  assumed  geometry. 


Figure  A-l.  Schematic  of  gas  volume. 


Assume  one-dimensional  heat  flow  in  the  x-direction,  and  further  assume  that  heat  enters  the 
volume  by  natural  convection  from  the  wall  temperatures.  In  the  first  instance,  neglect  any  mo- 
tion of  the  interfaces,  and  compute  the  static  temperature  distribution.  For  an  infinitesimal  slice  of 
the  volume  of  length  dx,  the  heat  balance  is  shown  in  figure  A-2. 


h-  dx-*l 


Figure  A-2.  Heat  balance  in  differential  element. 
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APPENDIX  A 


The  heat  balance  is 


- (qx  + -to  dX) 


+ q = 0 , 


where 


= heat  conducted  in  x-direction  = -kvr2  dT/dx  , 
q = heat  convected  into  gas  = H2wr  dx  (Th  - T)  , 
k = thermal  conductivity  (W/mK)  , 
r = radius  of  pipe  (m)  , 

H = convective  heat  transfer  coefficient  (W/K/m2) 


The  governing  differential  equations  for  the  hot  and  cold  columns  are 


d2T  A 2H  , 

^ + kF  (Th'T> 


d T A 2H  , 

£5  + kF  (Tc-T> 


= 0 , 


A new  temperature  parameter  is  defined  as 


6h  = T - Th  . 


e = t - t 

C C 


In  addition,  define  a length  parameter  for  the  cold  side. 


and  let 


A2  = 2H/kr 


The  governing  equations  thus  become 


dx2 


a20. 


o , 


d20 
c 

dy2 


X20  =0  . 

c 


The  solutions  of  equations  (A-4)  to  (A-9)  are 

0h<x)  = Bj  sinh  Ax  + B2  cosh  Ax  , 

0 (y)  = B,  sinh  Ay  + Bu  cosh  Ay  . 
c 1 

The  boundary  conditions  are 
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( A— 8 ) 


(A-9) 


(A-10) 

(A-ll) 


T(x=0)  = T , 
n 

T(X=V*h)  = Tc  ' 

Th(x=*h)  = Tc(X=*h)  ' 

dTh(X°*h)  = dTc(X=th) 
dx  dx 


In  terms  of  the  new  parameters. 
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Equations  (A-12)  and  (A-15)  make  Ba  = B<  = 0.  If  now  lc  is  assumed  equal  to  l },  (on  the  average), 
then  by  symmetry  equations  (A- 13)  and  (A- 14)  may  be  replaced  if 


T(*=*h)  = 


T + T 
h c 

2 


Thus, 


and 


so  that 


(t  - T,  \ sinh  Ax 
a , , Vc  h) 

®h  X 2 sinh  U 

h 


e (y) 
c 


(Th  - Tc)  Sinh  Xy 
2 sinh  A l 

c 


(A-16) 


(A- 17) 


For  the  analysis,  average  values  may  be  used  so  that 


T 


h 


o 


and 


T 

c 


0c<y)  dy 


I 
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APPENDIX  A 

When  equations  (A- 16)  and  (A- 17)  are  substituted  and  the  indicated  integrations  are  performed, 
the  results  are 


T - T 

T = T + ..  „ C . , „■  (cosh  AH.  - l\  (A- 18) 

h h 2A£h  sxnh  A \ h J 

and 

T - T / v 

T = T + ttt-j -f,  (cosh  AS.  - 1)  . (A-19) 

c c 2X1  sinh  XI  \ c / 

c c 

These  average  temperatures  are  used  in  the  numerical  model. 

A more  sophisticated  model  assumes  that  the  hot  side  liquid  is  asymptotically  cold.  In  this 
model,  there  are  three  temperature  zones.  The  formulation  is  the  same  for  equations  (A- 18)  and 
(A-19),  except  that  the  simplifying  assumption  of  symmetry  is  not  used.  The  solutions  only  are 
given  here: 


T (x)  = B,  e*LX  + T 
LI  c 


T (x)  = B,  ^ — sinh  Ax  + (b  + T - T \ cosh  A x + T 
G V 1 c hi  G 


1 A, 


h ' 


T (x)  = Bc  sinh  A_(H  + i - x)  + T 


3.  * (t  - T ) cosh  \ i - 1 
1 B \ h c/ L G h 


c h 


+ tanh 


0] 

(Vc)  v»] 


-L  1 / 

B = sinh  A £,  + cosh  XI.  + - — tanh  XI  (XTcosh  A JL  + A_  sinh  XI. 
Ag  Gh  GhAG  Gc\L  Gh  G 


h) 


B.  = 


3l(^  sinh  Acih  + cosh  XcZhj  + (tc  - Th)(cosh  A^  - l) 


sinh  A„Jt 
G c 


and  the  average  hot  and  cold  temperatures  are 


Bj  + Tc  - Th 

cosh  AS.  - ll  + — — sinh  A„H 

G h / A-H, 


G h 


= 'h  + ^Th  (C°Sh  Vh  ‘ X)  + ~ 

n 

. , . (cosh  11  - 1^ 

c A„H_  V G c ) 


G*h  ' 


T = T + 
c 


G c 
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where 


T = temperature  of  liquid  < 
L 

G = gas  , 

£ = length  of  cold  column  , 

c 

£ = length  of  hot  column  , 

h 

B = constant  of  integration 
i 


Even  more  sophistication  may  be  introduced  by  dropping  the  one-dimensional  restriction.  Drop- 
ping it  leads  to  two  governing  equations, 


»2»i  (\  * *4  ./vdi  SA  , „ , 

ITT  * \—r~)  „2  \ ~ 7 » » 


where 


T - T 


i T - T 
h c 


5 ■ -/('h  * 'o)  ■ 

R = liquid  resistance, 
f = r/R, 
i = L,  G, 


with  theso  boundary  conditions  on  the  liquid: 


0r  (5,1)  = o , 

Li 

e L (—•*>,  y ) = o , 


ay 


(5,o)  = o , 


eT  ( o,y ) = e (0,y) 

L vj 
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and  with  these  boundary  conditions  on  the  gas: 
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, ({ „ J1- 0 1 ' ‘yk  * *c) 

° S I0'  5 s *h/(V  * lc)  ■ 
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G 

3f 


(C,0)  = 0 


3 0 3 6 

TT  (0'v)  = TT  (0'y) 


eGd,y)  = o 


These  equations  may  be  solved  by  separation  of  variables  and  lead  to  two  eigenvalue  prob- 
lems resulting  in  a double  infinite  series  in  hyperbolic  functions  and  Bessel  functions.  This  tedious 
solution  is  left  for  the  reader  if  he  should  desire  to  pursue  it. 


The  original  differential  equations  (A-2)  and  (A-3)  are  altered  when  gaseous  motion  is  con- 
sidered. Equation  (A-3)  becomes 


d?T  2H 
dx2  Rk 


4?c  u 


/ \ p dT 

(T.  - T) =0  , 

V h / k dx 


(A-20) 


where  u = liquid  viscosity.  If  one  assumes  that  the  initial  static  temperature  distribution  is 
developed  prior  to  any  volume  change  (gas  velocity  onset),  equation  (A-20)  may  be  used  to  ex- 
amine the  relative  importance  of  the  convective  and  conductive  terms.  Assume  that  a temperature 
distribution  exists.  Assign  order  of  magnitude  values  to  the  derivatives: 


T - T 

dT  _ h c 


dx 


■>  T - T 
d2T  h c 


dx-1 


For  the  conductive  term  d2T/dx2  to  be  negligible  compared  with  the  convective  term  (*Pc„u/k) 
dT/dx,  P 


A 
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or 


Vc  ui 

-g  h 


» 1 


(A-21) 


For  air, 


y =1.2  kg/m3  , 

= 1.05  x 103  w-s/kgK,  specific  heat  at  constant  pressure, 

l , = 0.1  m , 

h 

k = 0.026  W/mK,  thermal  conductivity  . 


For  the  velocity,  assume  sinusoidal  displacement  of  ±0.02  m at  3 Hz,  so  that  the  maximum 
velocity  is  2nf(0.02)  = 0.38  m/s.  Therefore, 


yc  un 
E 


h 


k 


1841  » 1 . 


This  indicates  that  during  oscillation  there  is  no  significant  effect  of  conduction;  consequently,  if 
the  gas  volume  moves  back  and  forth  as  a slug,  the  temperature  is  governed  only  by 


d6 

dx 


+ 


yc  ur 

E 

2H 


e 


0 , 


which  implies  that 


/ fc  uR  \ 

8 = B6  exp  ^ x ) ' 


and  if  the  boundary  condition  is  that  0 = 0,  then  Be  = 0,  and  the  temperature  remains  at  its  initial 
constant  value. 

In  reality,  the  volume  expands  and  contracts,  and  there  is  conduction  with  a lag  time,  so  that 
the  initial  temperature  distribution  eventually  decays  to  a more  uniform  overall  temperature  due 
to  the  mixing  effect  of  the  convective  terms.  In  addition,  more  complex  boundary  conditions  exist. 
A thermal  lag  through  the  wall  or  the  device  supplies  heat,  and  there  is  conduction  axially  in  the 
walls.  A considerable  effort  would  have  to  be  expended  to  thoroughly  analyze  the  complete 
problem,  given  the  velocity  distribution.  If  we  add  now  the  complexity  of  coupling  the  transient 
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three-dimensional  energy  equation  and  the  Navier-Stokes  equations  for  the  gas  to  the  time- 
dependent  oscillating  liquid  equations,  the  problem  becomes  massive. 

A tentative  time  to  solve  such  a problem  is  1 to  2 man  years.  Such  a solution  would  be  highly 
desirable  if  any  great  effort  or  expenditure  were  ever  contemplated  for  the  Fluidyne. 
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This  appendix  presents  the  Fluidyne  pump  computer  program  written  in  FORTRAN  IV  for  a 
360/70  IBM  computer.  Five  first-order  differential  equations  are  integrated  by  DVERK,  an  IBM 
supplied  Runge-Kutta  scheme  that  uses  a fifth-  and  a sixth-order  algorithm  to  reduce  the  global 
error  to  that  specified  by  the  user.  It  chooses  the  time  increments  within  the  desired  steps  of 
printout,  which  are  specified  by  XEND.  The  program  has  five  nested  parameter  DO  loops  that 
compute  the  response  of  the  Fluidyne  while  the  following  parameters  are  swept: 

Output  length 
Equilibrium  liquid  height 
Equilibrium  gas  volume 
Output  area 

Maximum  hot  temperature 

All  units  are  specified  in  SI  units,  except  for  the  gas  constant,  which  is  in  English  customary 
units,  but  is  converted  in  the  program.  All  outputs  are  in  SI  units. 

A sample  printout  is  given  for  case  1 of  table  I in  the  main  body  of  the  report. 
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Cl /C*/iE  1 2.  1 7 


HI1 


FlILYN  FORTRAN  PI 


NAVAL  P 1ST  jRALLATE  SCHLOl 


C THE  PROGRAM  FLUDYN  SCLVfS  THE  GOVERNING  EQUATIONS  FOR  THE  FLUIUYNE  3 FLUCuClO 
C FJ HP  AS  FORMULATED  dY  T.  OR/EaIECkI  ON  27  SEP  1977  FLG03C20 

DIMENSION  Y (51  ,C t 29  ) ,R( 5 ,9  1 FLOC0030 

DIMENSION  V PR  I ME  I 5 1 FLU03090 

CGMMCN  F9, PS.R95 , RHOL IQ. G. A29,H95.P3,P1 ,K 3 1,  Hi  3,  A 12, A 95  FLU  3 ">050 

EXIEFNAL  fcn  flucjceg 

C INPOT  PARAMETERS:  CCLD  TEMPERATURE ,K ; TEMPERATURE  INCREMENT, K;  FLoJGC 70 

C cwulLld~ldM  ht IGHT  S , M ; E QU l L lbKIUM  GAS  VOLUM  E , M3  I A MB  1.E  N T PR E S SUR E . KP A ; F L GCOCdO 
C COLD  LEG  A RE  A , M 2 ; HC  T LEG  ARE  A , M2 ; CUT  PUT  LEG  AREA,M2;CCLC  LEG  L ENG TH , MFL OC GC50 
C HOT  LEG  LENGTH, MlOuTPjT  LEG  L ENGTh , M : L I QUI 0 OENS I T Y, KG/M 3 ; FLGCOIOO 

G GAS  GCNSTAM ,FT-LBE/L3M-F;PI; L I QUID  V I S COS  I T Y, KG/ M/S ; M AX IMUM  HOT  LEG  FLUCC110 
G TcMPGRATLRE,K;GRAvI  TATICNAL  GCN S T A NT , M2 / S ; I I ME  INCREMENT ,S.  FLU00120 

RE  AC (5, 53  I T 1 ,OEL  T,  ht , VE  ,P5,  A1  2,  A29,  A95,  HC ,HH, HO.RHCL IQ  EL  LOO  130 

REACif ,51  ) R,PI , VIS, 1 M AX , G , DT  F1000190 

FLU001 50 


P5  = P5*  1 000 • 

FLUC01 60 

R-R *2.201* 1 . E*fa. E 5*1  2000. * (.  0259*  *3) 

FLU00170 

DO  60  JHC-1, 10 

- -FLU00180 

hC  =F  LCAT ( JHC) 

FL  UC  C 1 50 

aRITFI6,2C0)  HO 

FLU00200 

200 

FCRMATI*  H0-*,F7.1 ) 

FL LOO  2 10 

DC  61  JHE=1 ,9 

FLUC0220 

F E = • 1 *F  LOA  T ( JHE 1 

FL  UO J 2 20 

RR 1 T E ( 6 , 201  ) HE 

F L GO  32  93 

201 

FCRMATI*  HE=*,F7.2) 

FL  UCU250 

DO  62  JVE-1,9 

F L UC  326.3 

Vfc=.C03Dl*FLCAT ( JVE) 

FLUCC270 

aR 1 TE( t,2C2)  VE 

FL  UO  3 2 60 

202 

FCRMATI'  VE= • , F S. 5 ) 

FL  LU 3250 

CG  62  J A9  5 = 1 * 3 

FLGCC3C0 

A95=.CC01*FLOAT( JA95) 

FLU  30  3 13 

WRITEI6.203I  A95 

F L CO  33  20 

203 

FCRMATI'  A95=',F9.5) 

FLUCJ220 

DO  69  jTMAX-1,6 

FLU 3 02  50 

T MAX =295. *1C.*F LCAT ( JTMAX) 

FLUCC250 

ak  1 16(6,209)  TMA  X 

FL  UC 3260 

209 

FCRMATI'  TMAX=',F7.1) 

FL  LC  3 3 70 

Na  =5 

FLUCC360 

N=  5 

FL  LOO  2 53 

C IMI  1AL  CONDITIONS 

FLGC3900 

x = c.c 

FLUC0910 

v i n = c.c 

FLU0J9  20 

Y 1 2 I =0  .0 

FLL03920 

>13)  =hE- 

FLUCD990 

Y( 9 )=F  = 

FLUO 3950 

Y ( 5 I =FE 

FLUCC960 

T 2 = T 1 

FLU009  70 

C= l 1 t .»VE/A24 

FLL00930 

F = ! '•  Sh|  C I -1  . 1/2  ./C/S  INHI  0 1 

FLUCC950 

y.  : '"i  ii  -o. j 

FLU00500 

Y f - . ' i 3 

FLUO 35 10 

i.  i = :. l 

FLU03520 

in---  1 

FLUO 3 533 

CvAf-=F5<VE/R/Tl 

FLUCCE90 

C C 

L'  PcT-II-N  : AtSI  STANCES  ANC  lNDj;TA\CcS 

FL  U0C5  50 

C A 

NO  o t S * A 1 S . 

FLUCD560 

MR  11  E I a , 25  ) 

FLUC0570 

A A 1 1 t.3t) 

FLU035E3 

LC  13  K-  l.K  1 

F L UO  35  50 

R3i=5.*Fl*VlS*YI3)/A13**2 

FLOCC&CO 

R29=  c.  *PI*VI  S*Y(  <,1  / A 29**2 

FLU336 13 

R9i  = 5.*F  I*V  1S*1  FiU-Ht*Yi  5 1 ) /A95*«2 

FLU0:b20 

hi  2 = thCL  1 L * Y I 3 ) / A 1 3 

FLU036  33 

F29  = RhCL ( ;*Y( 9) /A  29 

ELU0:693 

h95  = Rh':Ut*IhO-hE*Y  IS!  I/A,5 

FLUCCt  50 

START  CC' ROTATION  CF  TIME  VARYING  SEQUENCE  AN 
X E NT  =FLCAT(K)*OT 
VI  = VE/2.-1Y (3 1-hE )*A13 
V2-Vt/2.-IYI9)-hti*A29 
V = v I ♦ V 2 


integration  of  the  ce-flucgcft 

FLU  C 36  7 3 
FLLCG6ED 
FLOCCC  5 3 
FL GOO  7 53 
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FILE:  FLUDYN  FORTRAN  PI 


NAVAL  POSTGRaCUATE  SCHOOL 


TC=T1+! T2-T1)*F 
TH  = T2*( Tl-T2)*F 

P1=GMASS*R*TC*TH/(V1*TH*V2*TC  I 
P2  = P1 

IF  I V ( II .GT.0.01  GC  TC  6 

ei  = i .0 

GO  TO  7 
81=0.0 
COST 1NUE 

P4=H2A*!-YPRIME( 1)-YPRIME(2) )+P2fR2A«(-Y(1)-Y(  2)  l+P.FCLIO/ 
1AA5)**2*81»RH0LIU*Y( A)*G 
P3  = PA-KHCLli/.72*YI2 )*ABS(  Y(  2) )/ ( A2  A- AA5  ) **2 
CALL  DVERK(N,FCN,X,V,XEND,TOL, I NO , C , NW , W , I ER I 
IF ( INC.LT .O.OR. I ER.GT .0 ) GO  To  20 
P6=P1/10C0. 

POWER=Y(l)*(Y(l)/AA5)**2*RHOLIC/2. 

Y3=Y ( 3 ) -HE 
YA  = Y ( A I -HE 
Y5=YI5I-hE 
P7=P6-F5*.0C1 
VA=V-VE 

»RITE16«30)  Y(ll,POwER,Y3,YA,Y5,X,P7,T2,VA 
T2=T2+DELT 

IFIT2.LT.TMAX)  GO  TO  10 

T2=TMAX 

CONTINUE 

WR1TE16.S7)  TC.TH.R.C.F 
MR l TE (6  .88  I 

FORMAT!*  TC=*,F8.3,*  TH=',F8.3,*  R = * . F10. 3 . • 0=*,EIO.A,* 

lFOPM/T(*  ENC  CF  PROGRAM*) 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CCNTINUE 

STOP 

FORMAT! 12F6.C) 

FORMAT ! 6F6 .0 ) 

FORMAT! 13) 

FORMAT  !7E12. AfF7.lt  • *,E12.5) 

FORMAT ! * OUTPUT  FlCW  OUTPUT  PCWER  COLO  DISPL  HOT  OISPL 
1 SPL  TIME  GAS  PRESSURE  HCT  TEMP  GAS  VOLUME*) 

FORMAT! • 1M3/S ) I WATTS)  (M)  (M) 

1 (S)  ( KPAI  (K)  (M3)  •) 

CONTINUE 

MR  I T E ! 6 « AO  ) I ER i INC 
FOPMATl 21 3) 

STCP 

ENO 

SLeROUTINE  FCNIN.X.Y, YPRIME) 

DIMENSION  Y(3),YPRIME(5) 

CCMMCN  PA,P5,RA5,RHCLIQ.G, A2A.HA5 .P3.P1 .R31.H1 3, A12.AA5 
IF l !-Y(  1 ) -Y I 2) I.GT.G.C)  GO  To  11 

e2=i  .0 

GO  TO  12 
B2=C.C 

YPRIME! 1)  = IPA-P5-RA5*Y(1)-Y(5)*RH0L I G*G-RHUL 1 0/2 . * ( -Y < 1 )- 
lA2A/A2A*e2)/HA5 

YPRIME  !2)  = (P3-Pl-R31*Y(2)-Y(3)*RHOLIQ*G)/H13 

YPRIMEI3)=Y(2)/A13 

YPRIME!A)=I-Y(1)-Y(2))/A2A 

YPRIME(5)=Y(1 I/AA5 

RETURN 

ENC 


FLUCC710 
FLU00720 
FLUC0730 
FLUC07A0 
FLUO0750 
FLUU0760 
FLU00770 
FLU00780 
FLUCC750 
2.*IY(  D/FLUC0800 
FLU00810 
FLUCC820 
FLU00830 
FLU008AO 
FLUC0E50 
FLU00860 
FLUC0870 
FLU00680 
FLUC0890 
FLUC0900 
FLU00910 
FLU00920 
FLUC0930 
FLU009 AO 
FLUC0950 
FLUC0S60 
FLU00970 
FLUC0980 
F = * . E 10 .AFLUOC990 
FLL01000 
FLUC1C10 
FLUO 10  20 
FL  L01030 
FLUOIOAO 
FLUO  10  50 
FL  UC 1 CEO 
FLU01070 
FLU01080 
FL  UC 1 050 
FLUC11C0 
FL  UO 11 10 
OUTPUT  0 1 F LUC  1120 
FLU01130 
(M)  FLUC11A0 
FL  UO 1 1 50 
FLUO 1 160 
FLUC1170 
FL  UO 1 1 80 
FLU  C 1 1 50 
FLUC1200 
FLUO  1 2 10 
FLU01220 
FLUO 12  30 
F LUO  1 2 AO 
FLUC1250 
FLUG12fcO 
FLUO  1 2 70 
F L UC 1 280 
Y( 2) ) »*2/FLUC1250 
FLUC1300 
FL  UC  1 3 10 
FLUO 1320 
F L UO 13  30 
FLUC13A0 
FLUO  1 3 50 
FLLC13E0 
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